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We present a measurement of the branching fraction B{Df — » fi'^Vfj,) using a 548 fb~^ data sample 
collected by the Belle experiment at the KEKB e+e^ collider. The Ds momentum is determined by 
reconstruction of the system recoiling against DK^X in events of the type e^e~ — > D^DKX, Dl 
Ds"/, where X represents additional pions or photons from fragmentation. The full reconstruction 
method provides high resolution in the neutrino momentum and thus good background separation, 
equivalent to that reached by experiments at the tau-charm factories. We obtain the branching 
fraction B{Df — > ^i^v^) = (6.44 ± 0.76(stat) ± 0.57(syst)) • 10~^, implying a Da decay constant of 
fo, = (275 ± 16(stat) ± 12(syst)) MeV. 



PACS numbers: 13.20.-v, 13.20.Fc 

One of the important goals of particle physics is the 
precise measurement and understanding of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements, fundamen- 
tal parameters of the Standard Model (SM) . To interpret 
precise experimental results on decays of B mesons in 
terms of the CKM matrix elements, theoretical calcula- 
tions of form factors and decay constants (usually based 
on lattice gauge theory, see e.g. [H) are needed. Decays 
of charmed hadrons in turn enable tests of the predictions 
for analogous quantities in the charm sector. Measure- 
ments of charmed meson decay rates with an accuracy 
that matches the precision of theoretical calculations is 
thus necessary for checks and further tuning of theoreti- 
cal methods. 

The purely leptonic decay _D+ (the charge- 

conjugate mode is implied throughout this paper) is the- 
oretically rather clean; in the SM, the decay is mediated 
by a single virtual boson. The decay rate is given 
by 

r{Dt^i+.e) = ^fhmlmu^{l-^y\V.,\' , (1) 

oTT 771/ 

where Gp is the Fermi coupling constant, mi and niD, 
are the masses of the lepton and of the Ds meson, respec- 
tively. Vcs is the corresponding CKM matrix element, 
while all effects of the strong interaction are accounted 
for by the decay constant fn^- While the decay rate is 
tiny for electrons due to the strong helicity suppression 
and since the detection of r's involve additional neutri- 
nos, the muon mode is experimentally the cleanest and 
the most accessible one. Decays with electrons can be 
used to study the backgrounds. 

The analysis described in this paper uses data from the 



Belle experiment Q at the KEKB collider Q correspond- 
ing to 548 fb^^. We study the decay Df — > fi'^Vfj, using 
the full-reconstruction recoil method first established in 
the study of semileptonic D mesons 0] • Similar analyses 
have also been performed by the CLEO-c Q and BaBar 
experiments. 

The Belle detector is a large- solid- angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of time-of- flight scintillation counters (TOE), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to de- 
tect Kl mesons and to identify muons (KLM). The de- 
tector is described in detail elsewhere Q . Two inner de- 
tector configurations were used. A 2.0 cm beampipe and 
a 3-layer silicon vertex detector were used for the first 
sample of 156 fb~^, while a 1.5 cm beampipe, a 4-layer 
silicon detector and a small-cell inner drift chamber were 
used to record the remaining 392 fb~^ 

This analysis uses events of the type e"'"e~ — > 
D*D^''^K^''^X , where X can be any number of addi- 
tional pions from fragmentation, and up to one photon 
Q . The tag side consists of a D and a K meson (in any 
charge combination) while the signal side is a D* meson 
decaying to D^^. Reconstructing the tag side, and allow- 
ing any possible set of particles in X, the signal side is 
identified by reconstruction of the recoil mass, using the 
known beam momentum and four-momentum conserva- 
tion. 

Tracks are detected with the CDC and the SVD. They 
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are required to have at least one associated hit in the 
SVD and an impact parameter with respect to the inter- 
action point of less than 2 cm in the radial direction and 
less than 4 cm in the beam direction. Tracks are also re- 
quired to have momenta in the laboratory frame greater 
than 100 MeV/c. A likelihood ratio for a given track to be 
a kaon or pion, C{K,tt), is obtained by utilizing specific 
ionization energy loss measurements in the CDC, light 
yield measurements from the ACC, and time-of- flight in- 
formation from the TOF 0. We require C{K,Tr) > 0.5 
for kaon candidates. The momentum of the lepton candi- 
dates is required to be larger than 500 MeV/c. For elec- 
tron identification we use position, cluster energy, shower 
shape in the ECL, combined with track momentum and 
dE/dx measurements in the CDC and hits in the ACC. 
For muon identification, we extrapolate the CDC track 
to the KLM and compare the measured range and trans- 
verse deviation in the KLM with the expected values. 
Photons are required to have energies in the laboratory 
frame of at least 50 - 150 MeV, depending on the detect- 
ing part of the ECL. Neutral pion candidates are recon- 
structed using photon pairs with invariant mass within 
±10 McV/c^ of the nominal 7r° mass. Neutral kaon can- 
didates are reconstructed using charged pion pairs with 
invariant mass within ±30 McV/c^ of the nominal /v" 
mass. 

Charged and neutral tag-side D mesons are re- 
constructed in D ^ KniT decays with n = 1,2,3 (to- 
tal branching fraction ss 25%). Mass windows were 
optimized for each channel separately, and a mass- 
constrained vertex fit (requiring a confidence level greater 
than 0.1%) is apphed to the D meson to improve the 
momentum resolution. D*-candidates are not directly 
reconstructed: we construct the mass of the system re- 
coiling against DKX, using the known beam momentum, 
and require a value within ± 150 MeV/c^ of the nominal 
£>* mass [l^l- A recoil mass M-ceciX) defined as the 
magnitude of the four-momentum Pbcams — Py, for an 
arbitrary set of reconstructed particles Y . Pboams is the 
momentum of the initial e^e~ system. Since at this point 
in the reconstruction X can be any set of remaining pions 
and photons, there is usually a large number of combina- 
torial possibilities. It is reduced by requiring the presence 
of a photon that is consistent with the decay D* I?s7, 
where the Dg mass lies within ±150 MeV/c^ of its nom- 
inal mass [lo| . Further selection criteria are applied on 
the momenta of particles in the e^e~ rest frame; for the 
primary K meson the momentum should be smaller than 
2 GeV/c, for the D meson larger than 2 GeV/c and for 
the Dg meson larger than 3 GeV/c. The energy of the 
photon from D* Dgj in the lab frame is required to 
be larger than 150 MeV, irrespective of its polar angle. 
To further improve the recoil momentum resolution, in- 
verse 111 mass-constrained vertex fits are then performed 
for the D* and Dg, requiring a confidence level greater 
than 1%. After implying these selection criteria, the aver- 



age number of combinatorial reconstruction possibilities 
is approximately 2 per event. The sample is further di- 
vided into a right- (RS) and wrong-sign (WS) part. If 
the primary K meson is charged, both it and the D me- 
son are required to have opposite flavor (strangeness or 
charm respectively) to the D* , to be counted in the right- 
sign sample; all other combinations are wrong-sign. If the 
primary K meson is a Kg, the assignment is based on the 
relative flavor of the D and D* mesons alone. The flavor 
of the Z?* is fixed by the total charge of the X, assuming 
overall charge conservation for the event. 

Within this sample of tagged inclusive Dg decays 
(named D^-tags in the following), decays of the type 
Dg ^iv^ are selected by requiring another charged 
track that is identified as a muon and has the same 
charge as the Dg candidate. No additional charged par- 
ticles are allowed in the event. Remaining photons not 
used in the described reconstruction are allowed only if 
their total energy is smaller than l.O/m GeV, where m 
is the number of such particles. After these selections, 
in almost all cases only one combinatorial reconstruc- 
tion possibility remains. Figure [T] shows the mass spec- 
tra of Mrec (DKXj) (corresponding to the candidate Dg 
mass) and of MreciDKX^jj.) (corresponding to the neu- 
trino candidate mass). 

We define nx as the number of primary particles in 
the event, where primary means that the particle is not 
a daughter of any particle reconstructed in the event. 
The minimal value for nx is three corresponding to a 
e^e^ ^ D*DK event without any further particles from 
fragmentation. The upper limit for nx is determined by 
the reconstruction efficiency; Monte Carlo (MC) simula- 
tion shows that the number of reconstructed signal events 
is negligible for nx > 10. As the efficiency very sensi- 
tively depends on nx, it is crucial to use MC simulation 
that correctly reflects the nx distribution observed in the 
data. Unfortunately, the details of fragmentation pro- 
cesses are not very well understood, and standard MC 
events show notable differences compared to the data. 
Furthermore, the true (generated) value differs from 
the reconstructed , as particles can be lost or wrongly 
assigned. Thus the measured (reconstructed) n^ distri- 
bution has to be deconvolutcd so that the analysis can 
be done in bins of n'^ to avoid bias in the results. 

To extract the number of Z^s-tags as a function of 
in data, two dimensional simulated distributions in n^ 
(ranging from 3 to 8) and the recoil mass M,-cc{DKX^) 
are fitted to the RS and WS data distributions. The sig- 
nal shapes for different values of n^^ (ranging from 3 to 
9 [l3l) of the signal are modeled with generic MC simu- 
lation [l^l, which has been filtered at the generator level 
for events of the type e"''e" DgDKX. The weights 
of these components, wf" , z = 3, ...8, are free parameters 
in the fit to the data. As a model for the background in 
the RS sample, the WS data sample is used. The nor- 
malization constants between WS and RS (which vary 
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M?ec (DKX7^)/GeVV 



FIG. 1: Top: recoil mass spectrum for Ds-tags (for both RS 
and WS samples) . Bottom: spectrum of missing mass squared 
for Df fi'^Vfj, candidates for the selected data. Error bars 
represent the statistical errors. The dark-shaded areas show 
the fitted background, the light-shaded bands show the fit 
with systematic uncertainties. The vertical lines indicate the 
signal regions. 



with n^) are another six fit parameters. Since the WS 
sample contains some signal (« 10% of the RS signal), 
these signal components for different njr values are also 
included in the fit as independent parameters. As a cross- 
check, the fit has also been performed using generic MC 
RS-sample backgrounds, which gives a negligible change 
in the results. A further cross-check involved the ran- 
dom division of the MC sample into two halves, using 



the shapes of the first half to fit the signal in the second. 
The resulting weights as function of fit to a constant 
of 0.990 ± 0.046, which agrees well with the expectation 
of 1. The total number of reconstructed D^-tags in data 
is calculated as 



i=3 



(2) 



where N' 



MC,i 
Ds 



represents the total number of recon- 



structed filtered MC events that were generated with 



(regardless of the reconstructed n^) and 



the fitted weight of this component. 

To fit the number of fJ-Vfi events as a function 

of , two-dimensional histograms in and the recoil 
mass MrcciOKX^n) are used. The shape of the signal 
is modeled with signal MC distribution. As MC studies 
show, the background under the ^v^ signal peak consists 
primarily of non- I^s decays (« 18% of signal), leptonic 
T decays (where the r decays to a muon and two neu- 
trinos, w 7%) and semileptonic Dg decays (where the 
additional hadrons have low momenta and remain unde- 
tected, w 3.6%). Hadronic Dg decays (with one hadron 
misidentified as a muon) are a rather small background 
component (< 2% of signal). Except for hadronic decays, 
which are negligible, all backgrounds are common to the 
eVe mode, which is suppressed by a factor of 0(10^). 
Thus, the eVe sample provides a good model of the 
background that has to be corrected only for kinematical 
and efficiency differences. Including this corrected shape 
in the fit, the total number of fitted /ij/^ events in data 
is given by 



E 

i=3 



(3) 



where N'^^ 



represents the total number of recon- 
structed signal MC events that were generated in the 
i-th bin of (regardless of the reconstructed ) and 



is the fitted weight of this component. 



The numerical result for N'g"^ is 32100 ± 870(stat) ± 
1210(syst), that for N'^^ is 169 ± 16(stat) ± 8(syst). The 
statistical errors reflect the finite number of data signal 
candidates. The systematic errors are due to the limited 
statistics of WS data and MC signal and background 
samples. The errors were estimated by varying the bin 
contents of data and MC distributions and repeating the 
fits. By this procedure the non-ncgiigible correlations 
among the fitted weights were taken into account. 

As the branching fraction of — > /ij/^ used for the 
generation of MC events is known, the branching fraction 
in data can be determined using the following formula: 



^rcc 



Tyrcc 
^ */ii/ 



MC,rcc 



(4) 
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where BuciDs ^ fii^^,) 
number of reconstructed ^Vf^ 
weighted according to the fit to data, i.e 



0.51% and 7V^^''''''= is the 
events in MC simulation, 



N, 



MC,rcc 



MC,i 



(5) 



2 = 3 



The average efBciency for the reconstruction of Ds fJ-iy^ 
decays, e^t^, is not needed exp hcitly for the computation 
of the branching fraction The final result is: 

BiD, fiiy^) ■ 10^ = 6.44 ± 0.76(stat) ± 0.57(syst). (6) 

The quoted statistical error reflects the statistical un- 
certainty of the fltted weights wf'^ and , including 
their correlations. The systematic error combines the 
contributions due to the statistical uncertainties of data 
and MC background samples (0.29), the statistical uncer- 
tainty of the signal MC distribution (0.41), muon track- 
ing and identification efficiency (0.18) and possible differ- 
ences in relative rates of individual Dg decay modes be- 
tween MC simulation and data (0.19). Since the branch- 
ing fraction is determined relative to the number of Ds- 
tags, the systematic errors in the reconstruction of the 
tag side cancel. Differences in the neutrino peak reso- 
lution between data and simulation have been found to 
have a negligible effect on the systematic error. 

Figure[2](top) shows the branching fraction determined 
in bins of n^. The result is stable within errors in n^; 
note that the errors shown for the bins are correlated. 
As a cross check, also the branching fraction in a limited 
range < 6 has been determined as (6.54±0.76(stat) ± 
0.57(syst)) • 10"'^, which agrees well with the result given 
above. Figure [2] (bottom) shows our result in comparison 
with the PDG Il0| value and recent results from other 
experiments 0, [o|7 

In conclusion, we have studied events of the type 
e+e" ^ DlD^'^K^^^X,Dl with X = m:{-f) 

where the Dg is identified in the recoil of the remain- 
der of the event. Normalizing to this sample of _Ds-tags, 
the branching fraction of Dg — > fiVf^ was measured to be 
(6.44 ± 0.76(stat) ± 0.57(syst)) • lO'^, which is in good 
agreement with the current PDG value of (6.1 ±1.9)- 10^"^ 
[lO[ and also compatible with recent results from BaBar 
(6.74 ± 1.09) • 10-3 0] and CLEO-c (5.94 ± 0.73) • lO^^ 
[5|. Finally we obtain the decay constant , using Eqn. 
© (with \Vcs\= 0.9730 0) 

/d, = (275 ± 16(stat) ± 12(syst)) MeV. (7) 

A simple average of the decay constants following from 
the cited measurements has an uncertainty of around 
10 MeV. Recently an LQCD calculation of significantly 
improved precision was performed, with the result fo^ = 
(241 ± 3) MeV [3]. This value is somewhat lower than 
the experimental average and the comparison with the 
experimental results may point to some inconsistency be- 
tween the two. More precise measurements are needed 



o 

X 



Q 



o 

X 



Q 



12 

10 
8 
6 
4 
2 


9 
8 
7 



8 



PDG 



BaBar CLEO-c BELLE 



FIG. 2: Top: B{Ds t^Vp) as a function of rixi the final re- 
sult is shown as the dark-shaded region. For comparison, the 
PDG value and its error is shown as the light-shaded region 
in the background. Bottom: our result compared with PDG 
[Tol l and recent BaBar Q| and CLEO-c measurements not 
yet included in the PDG 2006 compilation; the dark-shaded 
region shows the weighted average of all measurements. 



for a firm comparison and will become possible in the 
near future at both B and tau-charm factories. 
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